Modulation Doping near Mott-Insulator Heterojunctions 
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We argue that interesting strongly correlated two-dimensional electron systems can be created by 
modulation doping near a heterojunction between Mott insulators. Because the dopant atoms are 
remote from the carrier system, the electronic system will be weakly disordered. We argue that the 
competition between different ordered states can be engineered by choosing appropriate values for 
the dopant density and the setback distance of the doping layer. In particular larger setback distances 
favor two-dimensional antiferromagnetism over ferromagnetism. We estimate some key properties of 
modulation-doped Mott insulator heterojunctions by combining insights from Hartree-Fock-Theory 
and Dynamical-Mean-Field- Theory descriptions and discuss potentially attractive material combi- 
nations. 
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I. INTRODUCTION 

The electronic properties of transition metal oxidesi*2*^ 
are determined by delicate balancing acts involving hy- 
bridized oxygen p and transition-metal d orbitals and 
strong correlations that suppress charge fluctuations on 
transition metal sites. Magnetic and transport proper- 
ties in these materials are extraordinarily sensitive to 
the character of the orbitals present at the Fermi en- 
ergy, and therefore to external influences like doping and 
strain. This sensitivity has motivated interest in the epi- 
taxial growth of oxide heterojunctions and artificial lay- 
ered oxides. This line of research seeks im part to emulate 
the achievements of semiconductor materials researchers 
who have over the past few decades learned to engineer 
the electronic properties of epitaxially grown semicon- 
ductor materials by exploiting lattice-matching strains 
and modulation doping. Because of the greater sensitiv- 
ity of their electronic properties and because of the wider 
range of phenomena (particularly magnetic phenomena) 
that occur, the implications for physics and for technol- 
ogy of substantial advances in the oxide case are likely 
to be enormous. In anticipation of future progress on 
the materials side, we explore in this paper some of the 
physics of modulation doping in epitaxially grown tran- 
sition metal oxides, emphasizing differences between the 
strongly-correlated-material case and the familiar semi- 
conductor heterojunction case. We find that strong cor- 
relations enhance the two-dimensional character of the 
metals that occur near modulation-doped heterojunc- 
tions and the range of doping which is possible without 
producing unwanted parallel conduction. We argue that 
two-dimensional electron systems produced in this way 
are likely to have remarkable properties and that modu- 
lation doping near interfaces between two-different Mott 
insulators may make it possible not only to create weakly- 
disordered low-dimensional strongly correlated electron 
systems, but also to engineer the compromises that oc- 
cur in these systems between different types of magnetic 
order. 

Transition metal oxides are prototypical strongly- 



correlated-electron systems. Understanding their elec- 
tronic properties has been one of the most challenging 
topics in condensed matter theory. Single-particle en- 
ergy scales like the widths of bands near the Fermi en- 
ergy are often comparable to or smaller than character- 
istic interaction energy scales, challenging band theory 
descriptions. Several different classes of transition metal 
oxides have been studied extensively revealing various 
interesting types of order involving spin, charge, and or- 
bital degree of freedomsi*^ and leading to fundamen- 
tal discoveries like high-T c superconductivity and colossal 
magnetoresistance. In the last decade, notable progress 
has been made in manipulating transition metal oxides 
by gating and by controlled layer by layer growth. Ahn et 
al. have applied the field-effect approach to ferroelectric 
oxide/high-T c cuprate heterostructures and successfully 
tuned superconducting properties^, Ohtomo et al& have 
observed unusual metallic behavior at Mott-insulator- 
band-insulator (MIBI) heterostructures realized by pre- 
cisely controlled growth of LaTiOa/SrTiOa layers. Very 
recently, Chakhalian 7 et al. have studied the interplay 
between magnetic and superconducting order at an inter- 
face between (La,Ca)Mn03 and YBCO. These achieve- 
ments not only provide new platforms for fundamental 
research, but also demonstrate the promise of devices 
with functionality that is based on the unique properties 
of strongly correlated oxide materials. 

There have also been important advances in the theo- 
retical description of artificially layered transition metal 
compounds. Efforts have been made2*2*i2*ii*i2i to un- 
derstand differences between surface and bulk proper- 
ties in strongly-correlated materials, providing insights 
into the main consequences of the absence of transla- 
tional invariance along certain directions. First principles 
calculations for PbTiOs/SrTiOs superlattices demon- 
strated interesting ferroelectric properties^ The ef- 
fect of spatially inhomogeneous multilayered structures 
on transport properties has been examined^ Recently 
Okamoto and Millisi^iSiii used combined insights from 
Hartree-Fock theory (HFT) and dynamical mean-field 
theory (DMFT) to investigate the LaTi0 3 /SrTi0 3 model 
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MIBI heterostructure systems mentioned above, and 
successfully described the interplay between long-range 
Coulomb interactions and strong short range correla- 
tions in the electron density distribution near MIBI het- 
erojunctions. They concluded that the unusual metal- 
lic behavior observed by Ohtomo et al. originates at 
the MIBI interface and that the properties of this in- 
terface are very different from those in the bulk, be- 
cause of an electronic surface reconstruction reminis- 
cent of the purely electronic^ reconstructions imposed 
by space-charge physics on systems with polar surface 
terminations^. 
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Motivated by this recent work we consider in this pa- 
per modulation doping near an interface between two- 
different Mott insulators, a MIMI heterojunction. The 
model system that we have in mind is sketched in Fig. 
[ |5|. Most classes of transition metal compounds are 
either ternary or quaternary, with additional specta- 
tor atoms that donate electrons to hybridized transition 
metal-oxygen orbitals near the Fermi energy. These sys- 
tems can be doped by replacing the spectator atoms by 
atoms with a different valence. Modulation doping of a 
MIMI heterojuntion is achieved by doping the larger gap 
material at a spectator atom location that is removed 
from the heterojucntion. The extra electrons then enter 
the upper Hubbard band of the lower gap Mott insulator, 
creating a two-dimensional doped Mott insulator that is 
trapped near the heterojunction by space charge elec- 
tric fields. The spatial separation between dopants and 
the carriers that reside in the upper or lower Hubbard 
bands should give rise to strongly correlated metals that 
are relatively free from disorder due to chemical doping, 
and are two-dimensional in character. These systems are 
illustrated schematically in Fig. [[I]. We study these sys- 
tems using both HFT and DMFT as in previous studies, 
and also demonstrate that a generalized Thomas-Fermi 
theory (TFT) can be employed to capture key qualita- 
tive physics of strongly-correlated heterostructures in a 
very direct way. TFT yields accurate results for charge- 
density profiles and for the critical doping 8 C associated 
with the onset of parallel conduction. We conclude that 
both the doping fraction 5d and the distance between 
the heterojunction and the doping layer play a role in 
the competition that occurs between different magneti- 
cally ordered states. 




FIG. 1: Modulation doping properties of different Mott- 
insulator Mott-insulator (MIMI) heterojunction classes. The 
figure illustrates the local electronic spectral function near 
the interface both before (left) and after (right) modulation 
doping. The upper Hubbard band spectral weight (yellow) is 
plotted with shading while the lower Hubbard band (red) is 
solid. Mott-Hubbard band bending near the interface is due 
to the electrostatic potential induced by the spatial separa- 
tion of dopants and carriers. The discontinuity in bands at 
the interface is determined by atomic scale physics particu- 
lar to an individual MIMI heterojunction. In analogy with 
semiconductor heterostructure terminology we define the fol- 
lowing classes of MIMI heterojunctions: (a) Type I: The Hub- 
bard gap of the smaller gap material is completely inside that 
of larger gap material. Both electrons and holes can then 
be trapped near the heterojunction, depending on doping. (b) 
Type II: The top of the lower Hubbard band or the bottom 
of the upper Hubbard band of the larger-gap material lies in- 
side the Hubbard gap of the smaller-gap material. Only one 
sign of carrier can be trapped near the heterojunction in this 
case, (c) Type III: The top of the lower Hubbard band or the 
bottom of the upper Hubbard band of the larger-gap material 
lies in the opposite band of the smaller-gap material. In this 
case charge transfer across the heterojunction occurs even in 
the absence of doping. Experimental determination of how a 
particular MIMI heterojunction system fits in this classifica- 
tion scheme is a key element of its characterization. In this 
article we study only Type I MIMI heterojunctions. 



In the next section we describe the single-band Hub- 
bard model used in this paper to address modulation- 
doped Mott-insulator heterojunction properties. In Sec- 
tion III we discuss results obtained for the electronic 
properties of this model using HFT, TFT, and DMFT. In 
section IV we discuss materials which might be suitable 
for modulation doping of Mott insulator heterojunctions. 
Finally in Section V we summarize our findings and spec- 
ulate on the potential of modulation doped Mott insula- 
tors. 



II. SINGLE-BAND HUBBARD MODEL 

The model system we focus on this paper is sketched in 
Fig. The heterostructure is composed of two different 
3d 1 perovskites, AMO3 and A'M'Oz, where both A and 
A' are group III elements and M and M' are group IV 
elements which have nominal 3d 1 electronic structure in 
this structure. Since the total number of electrons per 
unit cell is odd, compounds of this type must be Mott- 
type when they are insulators, unless translational sym- 
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FIG. 2: Schematic representation of the heterostructure stud- 
ied in this paper. We choose the x — y plane as the interface 
plane and z as the layer-by-layer growth direction. The sym- 
bol D (white circle) denotes a dopant layer containing a frac- 
tion 5d of dopant atoms with a valence larger than or smaller 
than the A (black circle) and A' (shaded circle) atoms. Our 
calculations are performed for a finite thickness film with N 
layers of transition metal M atoms (black dots) on each side 
of the heterostructure. In the one-band Hubbard model elec- 
trons hop between M sites only and are influenced by the 
space-charge field caused by the spatial separation between 
the dopant atoms and carriers in the upper Hubbard band. 



metries are broken. Modulation doping is achieved by 
replacing some of the A atoms in the larger gap insulator 
by elements with a different valence. In this paper we as- 
sume electron doping for convenience, although the hole 
doping case is completely equivalent, apart from (impor- 
tant!) materials specific details which we don't attempt 
to model in this qualitative study. If we assume that a 
fraction Sd of the A atoms in a single layer of the larger 
gap material is replaced by donor atoms then the sum 
over all M atom layers of the upper Hubbard band frac- 
tional occupancy must be 6d in order to accommodate 
the extra electrons. When modulation doping is success- 
ful the added electrons reside in the lower gap material, 
placing them some distance from the ionized donor atoms 
and reducing the importance of the chemical disorder 
normally associated with doping. 

For this qualitative study we use a single-band Hub- 
bard model that ignores any orbital degeneracies that 
might be present. Modulation-doping depends critically 
on the long-range Coulomb interactions so these must be 
realistically represented in the model. Our model Hamil- 
tonian includes hopping, short-range repulsion, and long- 
range Coulomb interaction terms H = H t + Hjj + HqouU 
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Zi = 1 for I = A, A' and (1 + S D ) for I = D. We 
do not account for randomness in the dopant layer in 
this paper. The index % denotes the position of a tran- 
sition metal ion (M) so that Ri = a(rii,mi,Zi) and 
Rf — a{rii + 1/2, + 1/2, Zi + 1/2) respectively in a 
perovskite unit cell with lattice constant a. For the sake 
of definitencss, we ignore the possibility of a d-band offset 
between the two materials, although these will certainly 
occur in practice. Given this assumption, a Type-I MIMI 
heterojunction will occur whenever the Hubbard U pa- 
rameter is large enough to produce insulating behavior in 
both materials. We consider a system with a finite num- 
ber 27V of layers labeled sequentially from left to right 
and define U(zi) — U\ for Zi = 1 to TV and U2 for z% = TV 
to 27V with U\ > U2 so that the larger gap material is on 
the left. We treat the Coulomb part of the interactions in 
a mean-field Hartree approximation. Since Coulomb po- 
tentials in the absence of doping are implicitly included 
in the model band Hamiltonian, in evaluating this poten- 
tial we include only the extra charges in the dopant layer 
and charges due to occupancy of lower or upper Hub- 
bard bands. To be specific, the reference background has 
charge per atom equal to —1 for each M site and +1 
for each A, A', and D site. As a result, the mean-field 
long-ranged Coulomb interaction is: 
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where pj = 'YliAtijcr) is the electron density on site j. 



III. TYPE-I MIMI HETEROJUNCTION 
ELECTRONIC STRUCTURE 

A. Hartree-Fock Theory 

In HFT the strong on-site Coulomb interactions is also 
treated in a mean-field approximation so that 



n^hii — > y^hj^) nig. 



(3) 



HFT is equivalent to minimizing the microscopic Hamil- 
ton in the space of Slater-determinant wavefunctions. As 
notedii previously there are typically a number of self- 
consistent solutions of the HF equations, corresponding 
to a number of local minima of the Hartree-Fock energy 
functional. The various minima usually are distinguished 
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by different types of magnetic order. Our philosophy 
in examining several different solutions without strong 
emphasis on their relative HF energies is that different 
types of order will occur near different interfaces but 
neither the single-band Hubbard model nor any of the 
electronic structure approximations we consider (or in- 
deed any known electronic structure approximation) is 
sufficiently reliable to confidently select between them. 
Indeed phase transitions between Mott insulator states 
with magnetic order and paramagnetic metallic phases, 
corresponding to magnetic and non-magnetic extrema of 
the Hartree-Fock energy functional, are often first order. 
(We will however make some conclusions of a more qual- 
itative nature concerning trends and tendencies related 
to modulation doping.) As explained more fully below, 
we find that the HFT electron density distribution near 
a MIMI heterostructure is sensitive mainly to the rela- 
tive orientations of electron spins on neighbouring metal 
sites on adjacent layers. Consequently, we present results 
only for usual bipartite antifcrromagnetic (AFM) and fer- 
romagnetic (FM) states, which in this respect cover the 
two possibilities. These two ordered states are metastable 
in both undoped and modulation doped regimes for the 
range of parameters we have studied. 

The results of our HFT calculations are summarized in 
Fig. We have chosen typical parameters for a one-band 
Hubbard model of perovskite transition metal oxides, 
taking Ui/t = 24, U 2 /t = 15, and U c = e 2 /eat = 
We can see from Fig. [2Jthat the modulation doping effect 
occurs for both AFM and FM states, although the details 
of the electron density distributions are quite different in 
the two cases. Short-range correlations therefore appear 
to play a relatively strong role in determining the charge 
distribution near MIMI heterostructures, in contrast to 
the MIBI heterojunction case in which they plaj*i^ a rel- 
atively minor role. The upper Hubbard band electrons 
are noticeably more confined to the interface in the AFM 
state case and spread further into the smaller-U layer in 
the FM state case. This difference in density-distribution 
follows from a corresponding difference in the compro- 
mise between band-energy minimization and interaction 
energy minimization in the two-cases. The ferromagnetic 
state which has all spins parallel maximizes the hopping 
amplitudes between sites, doing a better job of minimiz- 
ing band energy at a cost in interaction energy. The 
bandwidth of the mean-field quasiparticle states is ~ t 
for FM states and ~ t 2 /U for AFM states. Increased 
doping should favor FM states over AFM states, at least 
within HF theory. From a real-space point of view, dop- 
ing frustrates the staggered moment order of the AFM 
state more strongly than it frustrates the FM order be- 
cause of the nearest- neighbor hopping term H t . In other 
words, doping favors the FM state over the AFM state. 

The doped electrons have a strong tendency to accu- 
mulate nearly completely in one layer in the AFM state 
case. Larger setback distances for the dopant layer should 
result in larger space-charge fields at the heterojunction 
and less opportunity for electrons to spread out away 
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FIG. 3: Electron density distributions from HFT for (a) AFM 
and (b) FM states as a function of Sd- The parameters used 
are Ui/t = 24, U 2 /t = 15, U c = 0.8, and TV = 5. 2 is the layer 
index for M site so that U(z) = Ui for 2 = 1 — 5 and U2 for 
2 = 6 — 10. The dopant layer is at zd = 3.5. 



from the interface, robbing the ferromagnetic state of the 
extra stability that it gains from the third dimension. We 
expect therefore that for a given doping level Sd, antifer- 
romagnetism will be favored by a larger setback distance 
for the dopants. A larger setback distance also favors 
the development of parallel conduction channel. These 
trends can be seen in the ground-state phase diagram 
plotted in Fig. ^ In summary, modulation doping in 
MIMI heterostructures may make it possible not only to 
create weakly-disordered low-dimensional strongly corre- 
lated electron systems, but also to engineer the compro- 
mises that occur in these systems between different types 
of magnetic order. 

We note in Fig.[3]that for the FM state, a parallel con- 
duction channel starts to appear adjacent to the doping 
layer at Sd = 0.375. For the parameters we have chosen 
modulation doping successfully places the carriers in a 
more remote layer up to this doping level. 
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FIG. 4: Ground state (T = 0) HFT phase diagram vs. doping 
concentration Sd and setback distance (zn+i — zd). AFM- 
NP denotes the antiferromagnetic state without a parallel 
conduction channel, and FM-NP (FM-P) denotes the ferro- 
magnetic state without (with) a parallel conduction channel. 
Larger setback distances favor antiferromagnetism over ferro- 
magnetism and strengthen the tendency toward development 
a of parallel conduction channel. 
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The HFT results can be understood using a Hubbard- 
model version of Thomas-Fermi theorjiSi. The TF equa- 
tion for this system are: 

p(p(z)) + vh{z) = const (4) 

where p(p) is the chemical potential at density p without 
long-ranged Coulomb interaction and vh(z) is the elec- 
trostatic potential for zth layer obtained from the charge 
density by solving the Poisson equation. In principle, 
p,(p) should be obtained from the exact solution of the 
three-dimensional Hubbard model. This input is unfor- 
tunately still unavailable. Instead, we can use HFT to 
obtain p(p). In this way we have separate versions of 
the TF equations for AFM, FM, and PM states. As for 
vh{z), in the continuum limit each layer can be approxi- 
mated by a 2-d uniformly-charged plane so that we have: 

=Sd\z~ zd\ - £ (p(z') -i)\z- A (5) 

where zd is the layer index for of the dopant layer and z' 
is summed over all electronic layers. Fig. [S] shows results 
calculated using this TFT for the same parameters as 
used in Fig. [31 The total electron density distributions 
are almost identical to those obtained from the full mi- 
croscopic HFT. We do note that the parallel conduction 
channel in the FM state appears at lower doping in TF 
theory than in the microscopic theory. 



FIG. 5: Total electron density distributions from TF theory 
for (a) AFM and (b) FM states with the same parameters as 
used in Fig [3] The results are close to those in Fig. [3] except 
for differences in the critical doping S c at which modulation 
doping starts to fail. 



The local-density approximation for correlations im- 
plicit in the Thomas-Fermi theory is obviously least reli- 
able in judging the relative chemical potentials for adding 
carriers to spatially separate layers. The discrepancy also 
occurs partly because the long-ranged Coulomb interac- 
tion was evaluated using a three-dimensional lattice ver- 
sion in the microscopic HFT whereas a continuum ap- 
proximation for the layers was used in the TFT calcula- 
tions. We expect that the TFT model is too simplified to 
determine the critical doping for parallel conduction 5 C 
accurately for any particular heterojunction, but it may 
be used to analyze trends. 

The most important consequence of strong local in- 
teractions in the Hubbard model is the emergence of a 
jump in the chemical potential when the electron den- 
sity per site crosses from less than one to more than 
one. As illustrated in Fig. the opening of the Hub- 
bard gap is accompanied by slower dependence of the 
chemical potential on density just above and just be- 
low p(z) — 1, i.e. by an increase in the thermody- 
namic density of states within the Hubbard bands. To 
capture these features we approximate the chemical po- 
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tential in the upper Hubbard band near p(z) — 1, by 
H{p{z)) ~ E c {U(z)) + (p(z) - l)/V(U(z)), where E C (U) is 
the bottom of the upper Hubbard band, and T>(U) is the 
thermodynamic density of states averaged over the en- 
ergy range of interest near the bottom of the band. This 
notation is chosen to emphasize similarities to semicon- 
ductor heterojunction physics. Using this result in each 
layer we find that 



E c (Ui) - E C {U 2 ) 
4irU c (z N+1 - z D ) 



O 



(6) 



where zjv+i is the index of the first metal layer on the 
small U side of the heterojunction. This simple and ap- 
proximate expression emphasizes that 8 C increases with 
E c (Ui) — E C (U2), decreases with U c and, as in the semi- 
conductor case, decreases with the donor layer set-back 
distance. ^From HFT we estimate that E C (U) « U for 
AFM states while E C (U) &U-6t for FM states. In both 
cases E c (Ui) — E C (U2) is approximately Ux — U%. As illus- 
trated in Fig. El T> tends to be larger for AFM states than 
for FM states. The precise form of the thermodynamic 
density of states near the band edge in any particular ap- 
proximation can only be determined numerically. This 
simple expression does not fully capture the difference 
between AFM and FM states, but it does capture some 
simple but important properties. More effective mod- 
ulation doping will occur materials combinations with 
larger U difference, and smaller U c (i.e. larger dielectric 
constant e) values. Additionally, because of stronger ten- 
dency to confine electrons in one layer in AFM state, S c 
is larger in AFM than in FM states in general. These 
features are confirmed by our numerical calculations. 

TFT is successful because the dominating energy scales 
are the electrostatic energy and the correlation energy 
arising from local correlations. The ground state electron 
density distribution is a result of competition between 
these two energy scales, which is accurately captured by 
the TF approximation. 



C. Dynamical Mean-Field Theory 

HFT provides a particularly poor description of para- 
magnetic (PM) strongly correlated states because it is 
unable to capture correlated quantum fluctuations. In 
the limit of large U it is clear that the thermodynamic 
properties (fx(p) for example) of a paramagnetic state 
are much more similar to those of ordered states than 
suggested by Hartree-Fock theory. To obtain a bet- 
ter description of paramagnetic modulation-doped elec- 
tron systems we appeal to dynamical mean-field-theory 
DMFT—. Full DMFT calculations are however too time 
consuming, even for our relatively simple model, and we 
therefore employ the two-site method, a minimal realiza- 
tion of DMFT that is able to capture the essential physics 
near a Mott transition— The two-site method has been 
used previously to describe a MIBI heterostructurei^. 




FIG. 6: Chemical potential /i versus electron density p for 
AFM, FM and PM HFT states of the one band Hubbard 
model with U/t= 24. The jumps at p = 1 for AFM and FM 
states signal the opening of charge gap. The exact chemical 
potential is likely intermediate between the AFM and FM 
HFT values. 



Following the general framework of DMFT and the no- 
tation used by Okamoto et aZ— , the electron Green's 
function for each in-plane momentum can be written 
as: 



G(z,z',k\\;w) 



H, - H, 



eff 



Coul 



Y,(z,z',lo) 

J (7) 

where H(Joul 1S gi ven m Eq. [21 The self energy 
T,(z, z ,u>) = 6 ZtZ 'Ti(z,u>) is obtained by solving a two- 
site quantum impurity model for each layer and satis- 
fying a set of self-consistency equation o 15 ! 23 . We note 
that the two-site method predicts that the critical value 
of U for the metal-insulator transition of a 3-d single 
band Hubbard model 11 to be U c w 14.7. The U x and 
U2 values we have chosen are both larger than U c so 
that both perovskites are Mott insulators in the two-site 
method. Fig. [7\ compares the results from DMFT and 
HFT for paramagnetic states. These results demonstrate 
that modulation doping is possible without magnetic or- 
der in DMFT. The failure of HFT in this respect is a well 
understood consequence of the importance of on-site cor- 
relation effects for MIMI heterostructure properties, and 
of the failure of HFT to capture these correlations. 

In Fig. [5] we plot DMFT layer-dependent electronic 
spectral functions 



A(z,u) = -(1/tt) / [dA||/(27r)^] Im G{z, z, k\\\ u -M0+). 

(8) 

Only the layer closest to the interface on the smaller-U 
side (z = 6 in the figure) develops finite spectral weight 
near the Fermi surface upon doping. The appearance of 
a peak in the spectral function near the Fermi energy 
in layers close to the interface is reminiscent of the find- 
ings of Okamoto et al~ for a MIBI heterojunction, who 
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FIG. 7: Total electron density distributions for the PM state 
calculated by (a) DMFT and (b) HFT. The DMFT results 
exhibit a modulation doping effect while those of HFT do not, 
implying that modulation doping near a MIMI hetero junction 
does not occur without on-site correlations. Note that Figure 
(b) can also be reproduced accurately by TF equation with 
the HFT u(p) of the PM state. 



refer to this tendency as electronic surface reconstruc- 
tion. The robustness of this phenomenon beyond the 
two-site method is not certain at present, nevertheless it 
is intriguing that it occurs in two quite different hetero- 
junction systems. This finding has a natural interpreta- 
tion in DMFT. In Eq. /x - H e c s J ul {z) acts like "layer- 
resolved chemical potential," which determines the to- 
tal electron density in layer z. For each layer one must 
solve a separate quasi two-dimensional quantum impu- 
rity model whose solution shows insulating (metallic) be- 
havior for electron density close to (away from) 1. The 
self-consistency equations ensure that all solutions are 
related. Consequently, DMFT generally predicts that a 
layers with electron density away from 1 (layer 6 in the 
present calculation) is metallic. 
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FIG. 8: Local spectral functions for each layer calculated 
from DMFT for 8 D = 0.0625 (solid line) and 0.3125 (dotted 
line). Only the layer closest to the interface on the smaller-U 
side (z = 6) becomes metallic upon doping. 




FIG. 9: Total electron density distributions from TF the- 
ory for AFM, FM, and PM states for the MIBI heterostruc- 
ture with U/t — 16, and U c = 0.8. The three-unit-cell-wide 
crossover is remarkably reproduced. The similarity in electron 
density distribution between these different states demon- 
strates the relatively weak dependence of MIBI heterojunc- 
tion properties on details of the on-site correlation pointed 
out by Okamoto and Millis. The dependence of growth direc- 
tion density distribution on local correlations is stronger for 
MIMI heterostructures. 
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D. Thomas-Fermi Theory for MIBI 
Heterostructures 

Our TFT can accurately reproduce the electron den- 
sity distribution near a MIBI heterostructurei^ calcu- 
lated in previous work. Following the notation of 
Okamoto and Millis, the electrostatic potential in the TF 
relation of Eq. 0]can be written as: 

V ^ = 52\z-z a ,\-^2p(z')\z-z>\ (9) 

C Z A> Z'^Z 

where za> — ±0.5, ±1.5, • • • , ±4.5. Fig. shows solu- 
tions of the TF equations for AFM, FM, and PM states 
which are in accurate agreement with Okamoto's results. 
The three-unit-cell wide crossover is clearly seen for each 
state. The reason for the weaker dependence on the de- 
tails of on-site correlations is that many layers have elec- 
tron filling smaller than one. Electrons only fill in the 
lower Hubbard band. As a result, the correlation gap 
between the Hubbard bands, which is sensitive to the 
details of the on-site correlation, does not affect the elec- 
tronic structure very much. 



IV. MATERIALS CONSIDERATIONS 

It is still a challenge to determine model parame- 
ters appropriate for particular transition metal oxides, 
although a number of different approaches have been 
proposedi2i2£i2& Although not fully mutually consis- 
tent, these ideas do provide a general picture of how im- 
portant model parameters vary, and are a useful rough 
guide to possible materials combinations that might ex- 
hibit MIMI modulation doping. RMO3 materials (R: 
Rare earth) appear to be an attractive possibility be- 
cause of their relative simplicity. It has been shown that 
YMO3 has stronger electronic correlation than LaMC>3 
because of the smaller tolerance^ factor /. Therefore 
YMO3 /LaMC>3 heterostructure appear to be a good can- 
didate for realizing modulation doping. In particular 
LaTiC>3 and YTiC>3 are both Mott insulators with dis- 
torted perovskite structures (GdFeC>3 type) having gaps 
w 0.2 eV and ~ 1 eV respectively^. They might be used 
to realize a modulation doped heterostructure if YTiC>3 
could be doped. We emphasize that at present we do not 
know how the t 2g d-bands are aligned at a heterojunc- 
tion between these two materials. Indeed, aside from the 
band-lineup issue, it is importnat to recognize that the 
simple model considered in this paper is not sufficiently 
rich to capture all aspects of the interface physics that 
can be relavant to modulation doping and to magnetic 
order in the interface layer. For example orbital degener- 
acy plays a key role in the magnetic state of bulk YTiC>3 
and in all liklihood would also play a role in determining 
the magnetic state of any two-dimensional matallic layer 
at the interface. 



Although our one-band model is intended to qualita- 
tively describe 3d 1 systems with cubic perovskite struc- 
ture, some of our results should be generalizable. As 
emphasized above, the modulation doping effect is a con- 
sequence of Coulomb space-charge fields and on-site cor- 
relations. We therefore do not expect that perfect affinity 
to the ideal perovskite structure to be of key importance. 
Other RMO3 type Mott insulators with i?=rare earth or 
alkaline earth and M =Mn,Cr might also be good can- 
didates, although there will certainly be additional com- 
plications because of the larger c£-valence that are not 
addressed at all in this work. Building up more real- 
istic models for potential building block materials is an 
important challenge for theory. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have presented some theoretical con- 
siderations related to modulation doping near hetero- 
junctions between two different Mott insulators, combin- 
ing insights from Hartree-Fock Theory, Thomas-Fermi 
Theory, and Dynamical Mean-Field-Theory approaches. 
Using typical parameters within a simple single-band 
Hubbard model, we predict that modulation doping is 
possible with doping layers set back from the hetero- 
junction by several lattice constants. Modulation dop- 
ing can be used to create 2-d strongly correlated elec- 
tron system with weak disorder and controllable densi- 
ties. These systems could prove to be an interesting plat- 
form for systematic studies of strongly correlated system. 
Indeed, we find that the magnetic phase diagram can 
be altered not only by the dopant density Sd but inde- 
pendently by the dopant layer setback distance. Unlike 
the case of MIBI heterostructures in which the Coulomb 
field dominates, the growth direction electron density dis- 
tribution in MIMI heterojunction systems also depends 
strongly on the character of in-plane ordering. In the 
AFM state electrons are more confined to the interface 
while in the FM state electrons spread further away from 
the heterojunction into the smaller gap material. The 
results of HFT can be reproduced remarkably well by 
Thomas-Fermi theory, indicating that the electrostatic 
energy and the correlation energy resulting from local 
fluctuations dominate the physics of the heterostructure. 
From the TF equation, we estimate the critical doping S c 
at which modulation doping starts to fail. DMFT calcu- 
lations show that modulation doping can occur without 
magnetic order , and that it requires only the on-site 
correlations that lead to the Mott-Hubbard gap. Layer- 
dependent spectral functions calculated using DMFT in- 
dicate that only the interface layer is metallic, reminis- 
cent to the earlier findings of Okamoto et alm^. 

Doped Mott insulators typically appear to have exotic 
properties when the doping is small and more conven- 
tional properties when the doping is large and the total 
band filling is well away from one, the value at which lo- 
cal correlations have maximum importance. In the case 
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of the extremely heavily studied cuprate systems, for ex- 
ample, this crossover is interrupted, by high-temperature 
superconductivity. It is interesting to consider whether 
or not the two-dimensional electron systems considered 
in this paper are Fermi liquids. Whereas bulk doping 
often leads eventually to a first order transitions be- 
tween a doped Mott insulator and a relatively conven- 
tional metal, modulation doping in a single or several 
layers may make it possible to realize high-density, low- 
disorder, two-dimensional exotic metals which carry re- 
flect the heritage of the three-dimensional Mott insula- 
tors from which they emerge. Since the very existence 
of these two-dimensional electronic systems depends on 
gaps that are entirely due to electron-electron interac- 
tions, it is clear that they cannot be adiabatically con- 
nected to non-interacting electron states. On the other 
hand, in the HFT description the doped state is a Fermi 
liquid with well defined quasiparticles. This approxima- 
tion neglects quantum fluctuations of the magnetic state 
however, and its predictions for quasiparticle properties 
may not be reliable. 

We have also speculated briefly on materials combi- 
nations that might be attractive to realize the physics 



discussed in this paper. Predictions of the phase dia- 
gram for particular materials combinations will require 
much more detailed modeling, and may be assisted by in- 
sights from experiment as well as from ab initio electronic 
structurei2i2& calculations. A detailed description would 
require many realistic features of perovskite materials to 
be addressed, for instance lattice distortions^*^ which 
may vary with proximity to the interface, and related or- 
bital degeneracy issues. Progress will require progress in 
materials growth and characterization and interplay with 
ab initio and phcnomcnological modelling. 
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